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ABSTRACT 

It has recently been suggested that the power spectrum of redshifted 21 cm fluctuations 
could be used to measure the scale of baryonic acoustic oscillations (BAOs) during the 
reionisation era. The resulting measurements are potentially as precise as those offered 
by the next generation of galaxy redshift surveys at lower redshift. However unlike 
galaxy redshift surveys, which in the linear regime are subject to a scale independent 
galaxy bias, the growth of ionised regions during reionisation is thought to introduce 
a strongly scale dependent relationship between the 21cm and mass power spectra. 
We use a semi-numerical model for reionisation to assess the impact of ionised regions 
on the precision and accuracy with which the BAO scale could be measured using 
redshifted 21 cm observations. For a model in which reionisation is completed at z ~ 6, 
we find that the constraints on the BAO scale are not systematically biased at z > 6.5. 
In this scenario, and assuming the sensitivity attainable with a low-frequency array 
comprising 10 times the collecting area of the Murchison Widefield Array, the BAO 
scale could be measured to within 1.5 per cent in the range 6.5 < z < 7.5. 
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1 INTRODUCTION 

The imprint of baryonic acoustic oscillations (BAOs) on the 
mass power spectrum (PS) provides a cosmic yardstick that 
can be used to measure the dependence of both the angular 
diameter distance and Hubble parameter on redshift. The 
wavelength of a BAO is related to the size of the sound 
horizon at recombination. Its value depends on the Hub- 
ble constant, and on the dark matter and baryon densities. 
However, it does not depend on the amount or nature of 
dark energy. Thus measurements of the angular diameter 
distance and Hubble parameter can in turn be used to con- 
strain tli£j20SsJbl£_eTOluHan of^^ cosmic time 
(e.g. lEisenstein et al.lll99Sl ; lEisensteinll2002r i. 

Galaxy surveys have proved to be a powerful probe of 
the linear PS on large scales and therefore provide a means 
to measure the BAO scale (e. g. iBlake fc Glazebrookl | 2003| ; 
Seo fc Eisensteinl l200l 120051; iGlazebrook fc Blakej bood; 



Hu fc Haimanl 120031: ISeo fc Eisensteinl 120071: lAngulo et ail 



2008). Indeed, analyses of the galaxy PS from the SDSS 



and 2dFGRS surveys have uncovered a BAO signal (e.g . 
ICole et al. I l2005l ; lEisenstein et al. I l2005l ; iPercival et al"1l2007l ; 



lOkumura et a.1.1 120081 : iGaztanaga et alj|2008al lbh. providing 
incentive for deeper all-sky galaxy surveys, using telescopes 
like the SKA0, Pan-STARRS^ and the LSSlEI, to measure 
the BAO scale and hence the dark energy equation of state 
more accurately. 

Galaxy redshift surveys are best suited to studies of 
the dark energy equation of state at relatively late times 
(z < 3) due to the difficulty of obtaining accurate redshifts 
for a sufficiently large number of high redshift galaxies. Al- 
though detection of the Integrated Sachs- Wolfe effect puts 
some constraints on the integrated role of da rk energy above 
z ~ 1.5 [see, e.g., iGiannantonio et al.l (|2008l ) and references 
therein], we currently have very limited information about 
the nature of dark energy at high redshift. If dark energy 
behaves like a cosmological constant, then its effect on the 
Hubble expansion is only significant at z < 1 and becomes 
negligible at z > 2. In this case, studies of the BAO scale at 
low redshift would provide the most powerful constraints. 
However, as the origin of dark energy is not understood 
we cannot presume a priori which redshift range should be 
studied in order to provide optimal constraints on proposed 
models. Probes of dark energy at higher redshifts have been 
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suggested. These include measurements of the PS from a 
Lya forest survey which could potentially be used probe the 
evolution of dark energy through measurement of the BAO 
scale for redshifts as high as z ~ 4 (|McDonald fe Eisensteinl 
120071 ). Similarly, studying the temporal variation of high res- 
olution quasar spectra may probe the evolution of dark en- 
ergy through the Sanda ge-Loeb test in the window 2 < z < 5 
l|Corasaniti et al1l2007l ). Measurements of the BAO scale in 
the reionisation era (z > 6) are likely to be complimentary 
to the constraints at lower redshifts provided by these tech- 
niques. 

It has recently been suggested that observations of 
21 cm intensity fluctuations might pr ovide an additional av- 
enue to measure the mass PS (e.g. iMcQuinn et alj 120061 ; 
Bo wman et al 1 l2006l : iMao e i~aT1 l2008h . and that future low- 
frequency arrays cou ld be used to measure the BAO scale at 
a range of redshi fts (|Wvithe et al.ll2008l ; IChang et aLll2008l; 
IMao fe Wull2008l ). Redshifted 21 cm surveys are sensitive to 
neutral hydrogen regardless of whether it is part of a resolved 
object. As a result, at high redshifts such surveys may be 
more efficient in measuring the large-scale mass distribution 
than galaxy redshift surveys, which must identify a very 
large number of individual galaxies. Furthermore, spectro- 
scopic emission line surveys have the advantage of probing 
a precisely defined redshift interval, whic h is unlikely to be 
true of a traditional galaxy survey (see ISimpson fc Bridle] 
l2006h . 

Planned low-frequency arrays such as the Murchison 
Widefield ArrajQ (MWA) and LOFAP0are designed to mea- 
sure the PS of 21 cm fluctuations a t z > 6 in order to probe 
the reionisation era. IWvithe et al.l (|2008h have argued that 
while the first generation of low-frequency arrays will not 
have sufficient sensitivity to precisely determine the BAO 
scale, future extensions could have the sensiti vity to detect 
the B AO scale with promisingly small errors. IWvithe et al.l 
(2008) assumed a semi-analytic model which predicts a scale 
independent value for the effective bias between the lin- 
ear matter PS and the 21cm PS. However the 21cm bias 
is expected to be strongly scale dependent, especially to- 
ward the end of reionisation when many ionised bubbles 
percolate through the intergalactic medium (IGM), creat- 
ing an excess of power on scales larger than the charac- 
teristic bubble size (e.g. iFurlanetto et alj|2004 ). In this pa- 
per we explore whether the scale dependence of the 21 cm 
bias will compromise the ability of redshifted 21 cm o bser- 
vations to measure the BAO scale l|Mao fc Wul 120081 ). We 
begin by reviewing the expected BAO signal and 21 cm PS 
in Section We describe the observation of the 21cm PS, 
including a discussion of the noise considerations in Sec- 
tion [3] In Section [3] we describe our fitting procedure for 
recovering the BAO scale and present our error analysis, 
before concluding in Section [5] To calculate these errors we 
adopt a set of cosmological parameters similar to those de - 
rived from the third year WMAP data (|Spergel et al.ll2007h . 
namely Q m = 0.24, Q A = 0.76, h = 0.73, ag = 0.0.76 and 
n b = 0.042. 



2 THEORETICAL BACKGROUND 

The sound speed at recombination (s) determines the length 
scale of the acoustic peaks in both the cosmic microwave 
background (CMB) anisotropies and the linear matter PS. 
We refer to this length as the BAO scale. The co-moving 
horizon size at decoupling (r) may be determined from 
knowledge of £l m and Qt, and an und erstanding of th e 
physics governi n g a ba ryon-photon fluid l|Page et al.ll2003l ). 
iKomatsu et alj (|2008l ) calculate this length scale to be 
146 ± 1.8 Mpc for the cosmological parameters extracted 
from the fifth year WMAP data. The corresponding angular 
scale 9a ~ t/Da (zrcc) is closely related to the position of 
the first acoustic peak in the CMB PS, which occurs at the 
scale of the mode which has compr essed once at the time 
of decoupling (see IPage et al.l 120031 ) . Following decoupling 
the same acoustic scale is imprinted in the matter PS, al- 
beit with a different phase. This signature is expected to 
be preserved on large scales where non-linear gravitational 
effects are small. Therefore the BAO scale may be used as a 
standard ruler to test the geometry of the universe and the 
role of dark energy. 



2.1 Theoretical expectations for the redshifted 
21 cm power spectrum 

Neutral hydrogen radiates at 21 cm due to the hyperfine 
transition between the singlet and triplet ground states. As- 
suming a contrast between the kinetic temperature of the 
IGM Tfe and the CMB temperature Tqmb, and efficient cou- 
pling of the spin temperature T 3 to the 21 cm signal will 
mirror the underlying density of neutral hydrogen. The Lya 
and X-ray flux emitted by the first galaxies is expected to 
ensure these cond itions hold durin g most of the reionisa- 
tion era (see, e.g.. IFurlan etto 2006). A radio interferometer 
is sensitive to fluctuations in the brightness temperature of 
neutral gas. For gas of mean cosmic density the expected 
temperature brightness contrast in redshifted 21 cm emis- 
sion is 



5T b 



Tt — Tcm 
(1 + *) 



L (l- 
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where r is the optical depth of the neutral gas to 21 cm 
radiation and xm is the mass-weighted neutral fraction of 
hydrogen in the IGM. 

Allowing for fractional fluctuations in the baryonic mat- 
ter density S(x), ionised fraction 5 x (x), and peculiar gas ve- 
locity 5 v (x) — jjfe (x) (where is the peculiar velocity 
of the gasjf|, the spatial dependence of t he brightness tem - 
perature fluctuations may be written (e.g. IMao fc Wull2008r i 
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5T b {x) 



5f b [l - (1 - x m )(l + S x )](l + 5)(l - S v 

TCMB 



(2) 



Assuming the peculiar velocity effec t is described by the lin- 
ear theory result from iKaiserl (| 1987f ) . i.e. 5 v (k) = — f^ 2 S(k) 
(note that / — ► 1 in the high redshift limit) where n = k.n 
denotes the angle between the line of sight and the Fourier 
vector k, the a ngle dependent 2 1 cm power-spectrum may 
be written (e.g. IMao fc Wull2008l ) 



2a; H i(l - xm)Psx + 



P 2 i(k,z) = ST b 2 [(x 2 m P Sl 
(1 - xm) 2 Pxx) + 2f^i 2 (x 2 u Pss - Shi(1 - x m )Psx 

r2 4-2 7-, i 

+/ fJ- x m P S6 }. 



(3) 



Pss, Pxx and Ps x are the linear matter power-spectrum, ion- 
isation power-spectrum, and ionisation-density cross power- 
spectrum respectively. In general the 21 cm PS will dis- 
play complex angle and redshift dependent structure. The 
difference in the angular dependence modulating the cos- 
mological (Pss) and astrophysical (Ps x ,Pxx) components 
may be able to be exploited to con strain both cosmological 
and astrophysica l parameters (e.g. iBarkana fc Loebl [2005a; 
IMao fc Wul 120081 ). The spherically averaged 21cm PS may 
be more simply related to the underlying mass PS via 



P2i(k,z) = b 2 i(k,z) 2 P(k,z)D 2 (z), 



(4) 



where Pi is the primordial linear ma tter PS extrapolated 
to z = calculated using CMBfast (|Seliak fc Zaldarriaga! 
1996), and D(z) is the linear growth factor between redshift 
z and the present. In equation Q, the prefactor b 2 i(k,z) 
is a term which is analogous to galaxy bias (but which has 
units of mK). However in contrast to galaxy redshift sur- 
veys, the value of &21 can be scale dependent, even in the 
linear regime, as a result of the long-range, non-gravitational 
correlations introduced into the ionisation structure of the 
IGM by the formation of HII regions. Whils t pecu liar veloc- 
ity effects will be present. ISeo fc EisensteirJ (|2005r) find that 
the BAO signal is preserve d in t he matter PS a t z ~ 3 . Fur- 
thermore, ISeo et al.l (120081 ) and IShaw fc Lewis] (|2008l ) show 
that even the effect of non-linear redshift space distortions 
on the spherically averaged matter PS may be able to be 
corrected to high accuracy. 

Equation @ implicitly assumes that the gravitational 
dynamics governing the large-scale structures probed are 
well described by linearised fluid equations. Non-linear grav- 
itational dynamics will be increasingly important at low red- 
shift. However at the redshifts of interest here (z > 6) non- 
linear evolution is unlikely to play a significant ro le (see, e.g. 
ICrocce fc Scoccimarrol |200S| ; llVlao fc Wul 120081 1. IMao et all 
(2008) find that the cut-off in the maximum wave-number 
probed by redshifted 21 cm PS in the range 7 < z < 9 has 
little affect on the constraints on cosmological parameters. 



2.2 Model for the evolution of the 21 cm signal 

In the early stages of reionisation when most of the hydro- 
gen is still neutral, the bias between the linear matter PS 
and the 21cm PS, 621, is well approximated by a constant 
with units of temperature. However the bias is expected 



to depend on scale during the process of reionisation since 
galaxies preferentially form in overdense regions, resulting 
in early formation of ionised bubbles from an overdense 
IGM. This effect can b e modelled semi-analytically (e.g. 
IWvithe fc Moralesll2007h. or with full numerical calculations 
of re ionisation (e.g. iMcQuinn et al. I l2007al lbl; iTrac fc Cenl 
l2007fl . Recently, semi-numerical schemes have been devised 
which allo w efficient calculation of the large-scale ionisa- 
tion field (iMesinger fc Furlanettol 120071 ; IZahn et al.l 120071 ; 
iGeil fc WvitheT ^OOS). We use simulations co mputed using 
the method described in lGeil fc Wyithc (2008) to model &21 • 
We refer the reader to the discussion of the semi-numerical 
modeling in that paper, and to the discussion of the under- 
lyi ng semi-analytic model fo r the reionisation of the IGM 
in IWvithe fc Morales! (|2007h . Only a brief overview of the 
assumptions and characteristics of the model(s) is included 

here. 

IWvithe fc Morales! (|2007l ) describe a semi-analytic 
model for the density and scale dependent ionisation frac- 
tion of the IGM. This model assumes that a certain frac- 
tion of the mass in collapsed halos forms galaxies, and that 
these galaxies produce ionising photons with some efficiency. 
The ionisation due to galaxies is opposed by the density 
dependent recombination rate. An overdensity dependent 
rate of ionisation may then be calculated as a function of 
scale. For the examples shown in this paper, the model pa- 
ramete rs are tuned such that re ionisation is completed by 
z = 6. IWvithe fc Morales! (l2007h use this model to explore 
the 21 cm brightn ess temperature contr ast as a function of 
scale and redshift. IGeil fc Wvithel (|2008T) use these scale and 
overdensity dependent predictions for the ionised gas frac- 
tion of the IGM to generate realisations of the real-space 
density and ionisation fraction of a region of the IGM. In 
this paper we use the semi-numerical model to compute 
the 21cm bias 621- We note that our calculation does not 
include redshift space distortions l|Barkana fc Loebl [2005a') 
which are not captured by the model used in this paper. 
The free parameters include the star formation efficiency, 
escape fraction of ionising photons, and the minimum halo 
virial temperature for star formation in neutral and ionised 
regions of IGM (for which we take 10 4 K and 10 s K respec- 
tively). The semi- numerical simulation was performed in a 
cube with 256 3 resolution elements and of side length 3 Gpc 
(co- moving). 

In Figure [T]we show the model 21cm PS at z = 6.5, 7, 
7.5 and 8, corresponding to average volume weighted neutral 
fractions x m = 0.14,0.28,0.42,0.53. The left-hand panels 
show the ionised gas distribution from the semi-numerical 
simulations for a box with side length 3 Gpc (co- moving). 
The panels in the second column show the spherically aver- 
aged model 21cm PS (grey lines). The panels in the third 
column show the model 21 cm bias (grey lines). The far right- 
hand panels show the BAO component of the PS, defined as 
the difference between the 21 cm PS, and a reference 21 cm 
PS computed from a mass PS without baryons, Pi,„b(k). 
The model BAO component, b 2 21 (k)D 2 (z) [Pi (k) - Pi, nb (k)], 
is shown in grey. The model PS is reproduced in Figure 2. 



2.3 Fitting formula for 21cm power spectra 

Following ISeo fc EisensteirJ l|2005h we parametrize the fit to 
the semi-numerical model of the spherically averaged 21 cm 
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PS as 



P 2 i, fit (M =B 2 (k) D 2 (z) Pi(k/a) + AP(k), 



(5) 



where AP(k) = co + cik + 02k 2 is a term to describe anoma- 
lous power and B(k) = 6 + &i + ^2^ 2 + &3^ 3 is the fitted 
bias in units of mK. The "dilation" or "phase" parameter a 
allows for a dilation/contraction of the matter PS used to 
fit the 21 cm PS. In our case the input value of a is unity 
since we are assu ming we know the cosm ological parame- 
ters governing Pi . ISeo fc Eisensteinl l|2005h note that for the 
parametrization the PS in equation ([5]), the errors in the 
known values of the wave-number dilation parameter a di- 
rectly translate into the errors in the measured BAO scales. 
In Section \4. 21 we calculate the expected redshift dependent 
errors Aq for an observational campaign which is described 
in Section [3] below. 

The anomalous power term allows for the distortion of 
the m atter power-spect r um d ue to non-linear mode cou- 
pling. ISeo fc Eisensteinl <|2005l ) find that by subtracting a 
smooth function (AP) from the matter PS derived from 
their N-body simulations, they are able to recover the shape 
of the matter PS and thus the BAO signal. Our semi- 
numerical model is intrinsically linear, and so does not in- 
clude anomalous power. However we explore the possible 
effect of anomalous power on the measurement of the BAO 
scale by including the correction term AP in the fitted PS. 
Inclusion of this term approximates the effect of degeneracy 
between anomalous power and scale dependent bias. In this 
paper we present results using fits both with and without the 
AP term. Note that the degeneracy between the parameters 
governing the anomalous power and bias terms means that 
the polynomial fits to these functions cannot be interpreted 
physically when the AP term is included (see Section U). 
However we find that this degeneracy does not significantly 
hinder the determination of a. 

More detailed cosmological constraints may be ob- 
tainable if the line-of-sight and tra nsverse BAO scales 
could be measured independently [see Wyithe et al.l (12008 



and references therein, and also Okumura et al 



IGaztanaga et al.l (|2008al ) and IGaztanaga et al.l ( 2008bh for 



recent constraints from galaxy surveys]. Indeed, a low- 
frequency array is naturally suited to measuring both the 
line-of-sight and transverse components of the 21 cm PS (see 
Section [3]). Such an analysis would require modelling of the 
peculiar velocity effects which give rise to the angular depen- 
dence of the 21cm PS in equation ©. As mentioned, the 
model adopted for the 21 cm PS does not include peculiar ve- 
locities. However we note that above fitting formula could be 
readily generalised to characterise the ionisation and cross 
density-ionisation power-spectra needed to describe the an- 
gular dependence in equation (3). This would allow con- 
straints on the dilation parameter, and therefore the BAO 
scale, both perpendicular and parallel to the line-of-sight to 
be extracted from the 21 cm power-spectrum . The signal- 
to-noise ratio estimates in I Wyithe et al.l ()2008h suggest that 
the MWA5000 will be similarly sensitive to each of these 
scales. 



3 MEASUREMENT OF 21 CM POWER 
SPECTRA 

In this section we discuss the sensitivity of a low-frequency 
array to the 21cm PS. The sensitivity to the PS depends 
on the array design (e.g. total collecting area and antenna 
distribution) as well as the observation strategy (e.g. to- 
tal integration time). We consider a low- frequency array 
which would be a 10-fold extension of the MWA. We re - 
fer to this array as the MWA5000 (|McQuinn et all 120061 ). 
The MWA5000 would consist of 5000 antennae (with each 
antenna composed of a tile of 4 x 4 cross-dipoles) with effec- 
tive collecting area A e ~ 16A 2 /4m 2 , where A = 0.21(l+z) m 
is the observed wavelength. The dipoles will be sensitive to 
a wide range of frequencies which translate into redshifts 
spanning the expected extent of reionisation. However, we 
assume that the receiver will limit the bandpass to a selected 
Btot = 32 MHz window per observation. We assume that the 
antennae are distributed with constant density in a core of 
radius 80 m (within which the u-v coverage is very high) , 
with the remainder distributed with a radially symmetric 
density profile p oc r~ 2 out to a radius of 1 km (so that the 

maximum baseline is ~ 2 km) . 

We use the prescription from lBowman et al. I (|2006h for 
the 3d sensitivity of an array to the 21cm PS. The array 
sensitivity is expressed in terms of the wave-vector compo- 
nents which are orthogonal k± and parallel kit to the line- 
of-sight. The modulus of the wave-number is then given by 
k = |fc| = y/fcy + k\. The sensitivity at large k depends 

crucially on the angular resolution of the array (and there- 
fore the maximum baseline), the bandpass over which the 
measurement is made Av, the integration time tint, and the 
density of baselines which measure a particular fc-mode. The 
noise may be divided into two distinct components. The 
thermal noise component AP 2 i,jv is proportional to the sky 

' l + z > 2 - 6 1 



temperature, where T s ky ~ 2501 
of interest, and may be written 



K at the frequencies 



AP 2 i,jv(fe) = 



sky D'AD f V 
Avt in t n(k ± ) \A e 



(6) 



where D is the co-moving distance to the centre of the survey 
volume which has a co-moving depth AD. Here n(k±) is the 
density of baselines which observe a wave- vector with trans- 
verse component k±. In the denominator N c denotes the 
number of modes observed in a fc-space volume d 3 k. In terms 
of the fc-vector components, N c = 2irk±Ak±Ak\\ V/(2ir) 3 
where V = DAD(X 2 /A e ) is the observed volume. Note that 
this calculation of N c incorporates the fact that fc-modes 
which do not fit in the bandpass Av are not observable. 
The sample variance component of the noise is given by 
AP2i,sv(k) = P2i(k)/y/N c . The sensitivity will improve 
with the square root of the number of independent fields 
observed. The total noise may therefore be expressed by 



AP 2 l(fc_L,fc|| 



AP 2 i,iv(fcx, fey , z) + AP2i,sv{kx,k\\,z) 
V^ficidsBtot/A^ 



(7) 

For the purposes of constraining the parameters in equa- 
tion ([5]), we average equation (0 over the viewing angle to 
obtain the spherically averaged noise AP 2 i(fc,z). Note that 
this calculation assumes perfe ct foreground removal over a 
bandpass with width Av [see IMcQuinn etafl (|2006l ) for a 
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Figure 1. Model power spectra and fits assuming a third-order polynomial bias and no anomalous power. The left-hand panels show the 
ionisation maps corresponding to redshifts of z = 6.5, 7.0, 7.5 and 8.0. The panels in the second column display the corresponding model 
(grey) and fitted (black) 21cm PS, and the linear matter PS (thick grey line). The panels in the third column show the fitted (black) 
and model (grey) scale dependent bias. In the far right-hand panels we plot the difference between the model 21cm PS with and without 
baryons (grey), and the recovered 21cm PS with and without baryons (black), as described in the text. The grey error bars correspond 
to spherically averaged noise as described in Section [3] In each panel the vertical lines indicate the range of k values used for the fit. 



discussion of the relationship between foreground removal 
and bandpass width] . In Section [4] we explore the effect of 
the choice of Av on the errors in the measurement of a. 

We assume that the spherically averaged PS is mea- 
sured in rik = 32 logarithmically spaced intervals with end 
points in the range 0.05 Mpc -1 < k < 0.2 Mpc -1 . Wave- 
numbers below this range would not fit within the nominal 



bandwidth of ~ 8 MHz at the redshifts of interest. Above 
0.2 Mpc -1 the BAO signal becomes very weak and the mat- 
ter PS may be non-linear. The binning in fc-space, which 
is approximately equivalent to bin widths of Afc ~ fc/20, is 
motivated by the scale of the expected BAO features. For 
clarity the grey error bars in Figures 1 and 2 represent the 
uncertainty on the spherically averaged PS for larger spec- 
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Figure 2. Model power spectra and fits assuming a third-order polynomial bias and quadratic anomalous power for the same simulation 
data as Figure 1. The first column of panels display the model {grey) and fitted {black) 21cm PS, and the linear matter PS {thick grey 
line) at redshifts of z = 6.5, 7.0, 7.5 and 8.0. The panels in the second column show the fitted {black) and model {grey) scale dependent 
bias. The anomalous power component of the fit is plotted in the third column. In the far right-hand panels we plot the difference between 
the model PS with and without baryons {grey), and the recovered PS with and without baryons {black), as described in the text. The 
grey error bars correspond to spherically averaged noise as described in Section [3] In each panel the vertical lines indicate the range of 
k values used for the fit. 

tral bins with width Ak = fc/10. We assume that three fields 4 FITS TO 21 CM POWER SPECTRA 

are observed (iVfi e ids = 3) for 1000 hrs each. Given the field 
of view at these redshifts, this would correspond to survey- 
ing a few per cent of the sky. 

In this section we describe the results of fits to the modelled 
21cm PS using the parametrized form in equation ([5]). 
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4.1 \ 2 minimization 

We adopt the x 2 function as an indicator for the "goodness 
of fit" of our model, 



x 2 w = E 



P 2 l(k,z) - P21,flt(p, k,z) 



(8) 



where <7j = AP2i(k, z) and p = {Pj} = 
(a, bo, bi, bi, &3, Co, ci, C2) is a vector containing the model 
parameters. The index i refers to the discrete values of fe 
corresponding to each data bin (see Section [3J). We restrict 
our fitting to the range 0.05 Mpc^ 1 < k < 0.2 Mpc -1 . 

We hav e ada pted the Numerical Recipes routines 
l| Press et al.l [l992l which implement the Levenberg- 
Marquardt technique to iteratively find the parameter set 
which minimizes the x 2 function for a model which is non- 
linear in its parameters p. In Figures 1 and 2 we show the 
resulting fits for the cases with a third order polynomial 
bias (black lines), with and without a quadratic anomalous 
power term respectively. The panels in the second column of 
Figure 1 (and the first column in Figure 2) show the fitted 
spherically averaged 21 cm PS as well as the linear matter 
PS. The panels in the third column of Figure 1 (and second 
column of Figure 2) show the fitted 21 cm bias. In Figure 2 
the panels in the third column show the fitted AP term 
(note that the input term is zero since non-linear effects are 
not incorporated into the semi-numeric code we have used 
to generate the ionisation maps). The far right-hand panels 
in each of Figures 1 and 2 show the BAO component of the 
fitted PS, B 2 (k)D 2 (z)P l (k) + AP{k) -b 2 21 {k)D 2 (z)P Unb (k). 
The noise in the fitted BAO component is due to the noise 
in the model calculation of the bias 621. 

The fit to the 21 cm PS is very good in all cases. When 
the anomalous power term is included there is a trade-off 
between the power in the 21 cm PS associated with the bias 
B and that due to AP, indicating some degeneracy. However 
as we show in the next subsection, this degeneracy does not 
significantly bias the constraints on the BAO scale. It may 
of course be possible to alleviate the degeneracies between 
the bias and the AP parameters by putting some physi- 
cally/numerically motivated priors on the values of co,ci 
and C2. 



4.2 Uncertainty in the measured BAO scale 

After marginalising over the remaining fit parameters, we 
may interpret the error in a for our fits as the frac- 
tional uncertainty in t he measurement of the BAO scale 
Sco & Eiscnstc inl I2005T) . The Fisher matrix F is given by 



Tegmark et aill 19971 ) 



\ -> \ -> 3P21, fit 21, fit 

dp—' 



(9) 



where n p is the dimension of p and all derivatives are eval- 
uated using the best-fitting p. Since we expect the noise in 
21 cm PS to be Gaussian, we may invoke the Cramer-Rao 
theorem to estimate the la (68.3 per cent) confidence in- 
terval around the best-fitting value of a from F. Assuming 
that all other parameters are marginalised over this implies 
errors around the best-fitting values for a of 



(Pfl 1 ) 



-1/2 



(10) 



In Figure [3] we plot the the best- fitting values of a and 
the associated la error as calculated in equation (|10p for 
several redshifts in the range 6.4 ^ z 8 corresponding to 
average neutral fractions of 0.1 < xhi < 0.5 for our reion- 
isation model. The left panel of Figure [3] assumes a fitted 
21 cm PS with no anomalous power term, whereas the fit- 
ted PS used to generate the plot in the right panel includes 
the AP term. In each panel the grey horizontal line indi- 
cates the input value a = 1. Note that the recovered val- 
ues of a deviate from the input values at xm ~ 0.1 — 0.2 
[z ~ 6.4 — 6.7). However for the observing strategy and k- 
space binning chosen, this deviation is within the estimated 
la statistical error (although only just at z ~ 6.4 where the 
neutral fraction is ~ 10 per cent). As reionisation nears com- 
pletion (xm — > 0.1), the 21cm bias becomes most strongly 
scale dependent in the range 0.05 Mpc" 1 < k < 0.2 Mpc -1 . 
The increasing systematic error toward z ~ 6.4 is therefore 
a result of the (relatively) low noise level in the 21 cm PS at 
these redshifts combined with the increasingly strong evolu- 
tion of the bias function within the fitted range of k- values. 

We find that the statistical errors are slightly larger 
when we include the extra parameters required to describe 
the anomalous power (since the matrix F is not diagonal). 
We note that our estimate of the attainable errors on the 
phase/dilation parameter a are a minimum in the sense 
that (i) we have assumed perfect foreground removal across 
each bandpass with width 8 MHz for a total bandwidth of 
32 MHz, and (ii) we have invoked the Cramer-Rao theorem 
regarding the minimum la error on a parameter from the 
Fisher matrix of the best-fitting model. 

The bandpass over which foreground removal can be 
performed may affect the errors on the model parameters, 
even for a fixed total bandpass £?tot, because the removal 
of foreground on larger scales allows for more of the strong 
large-scale (small k) BAO peaks to be included in the fit. 
However we find that at the high redshifts that we are con- 
sidering, increasing Av above 8 MHz has little affect on the 
errors in the recovered values of a. On the other hand, we 
find that decreasing Av increases the systematic uncertainty 
in the recovered values of a toward the end of reionisation. 
For model fits without the AP term, the systematic error 
exceeds the statistical error when iehi < 0.2 (corresponding 
to redshifts z < 6.7), indicating that our maximum likeli- 
hood indicator is biased. This highlights the potentially sig- 
nificant advantages of including the larger scales accessible 
with wider bandpasses B. 

The systematic component of the errors will in general 
be sensitive to the functional form chosen to describe the 
21cm PS. We find that the effect of increasing the order 
of the polynomial used to fit the bias is somewhat degener- 
ate with the chosen value of n^. For example, with coarser 
k-space binning of nt = 16 the errors in a decrease by 
~ 0.5 per cent when the bias is modelled by a fourth or- 
der polynomial rather than a third order polynomial. This 
effect is not observed for nt = 32 (see Figures [3] and [4}. 
This implies that sampling of the BAO features must be 
sufficiently fine. However, the affect of nt on the systematic 
errors in a is only significant for xm < 0.2 when the bias is 
evolving strongly over the range of fitted k- values. 
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Figure 3. Fitted values and uncertainty for the a using a fit with third-order polynomial bias with {left panel) and without {right panel) 
a quadratic anomalous power term. The error bars represent the lcr (68.3 per cent) confidence intervals calculated using from the Fisher 
matrix by marginalising over all of the other parameters in the fit. 





Figure 4. Fitted values and uncertainty for the a using a fit with fourth-order polynomial bias with {left panel) and without {right 
panel) a quadratic anomalous power term. The error bars represent the 1<t (68.3 per cent) confidence intervals calculated using from the 
Fisher matrix by marginalising over all of the other parameters in the fit. 



5 SUMMARY AND CONCLUSIONS 

There has been recent interest in the utility of redshifted 
21cm fluctuations from the IGM a s a probe of the mass 
PS at a ranee of redshifts (e.g. [McQuinn et al.l 120061 ; 
Bowman et"alTl2006l: iMao et all l200Sl : IWvithe et all 120081 : 



Chang et a l. 2008; Mao fc Wul 120081 ). During the reionisa- 



tion era (z > 6), the appearance of HII regions generates 
scale dependence of the bias relating the 21 cm PS to the 
linear PS. In this paper we have set out to test the extent to 
which the scale dependence of the 21 cm bias will interfere 
with the ability to extract the BAO scale from the PS of 
21 cm fluctuations. 

We have utilised the results of semi-numerical models 
for the evolution of the ionised hydrogen distribution on 
large scales in order to construct an estimate for the 21 cm 
PS at high redshift. For a model in which the IGM is fully 
reionised by z = 6, and assuming the design parameters for 
a 10-fold expansion the the MWA, we find that the scale 
dependent 21cm bias limits the precision, but not the ac- 
curacy, of measurements of the BAO scale in the window 
6.5 < z < 8.0. We parametrize the bias and anomalous 
power terms in the 21 cm PS as polynomials, and for the tele- 
scope design assumed, find errors on the spherically averaged 
BAO scale of < 1.5 per cent in the window 6.5 ^ z ^ 7.5. At 
redshifts closer to the end of reionisation, the strong scale 
dependence of the 21 cm bias introduces systematic error in 



the BAO scale, which exceeds the level of the statistical un- 
certainty for narrower bandpasses {B < 8 MHz). We find 
that the sensitivity of a low-frequency array to the spheri- 
cally averaged BAO scale decreases with increasing redshift, 
becoming > 2 per cent by z ~ 8. 

The details of the quantitative results presented in this 
paper are sensitive to the reionisation model assumed. More 
generally, the constraints available at a particular redshift 
will be sensitive to the details of the reionisation history. 
The best constraints are likely to be derived from observa- 
tions of the Universe when the IGM is ~ 50 per c ent ionised 
- at w hich time the fluctuations are maximised dLidz et"ai1 
2008) - but the bias is relatively weakly dependent on scale 
at the scale of the BAO signal (and hence does not intro- 
duce a large systematic error). As a result, if reionisation 
completes earlier than z ~ 6, the most precise constraints 
will be obtained at higher redshift than reported here. How- 
ever, as the foregrounds become brighter when observing at 
higher redshift, an early reionisation scenario would decrease 
the maximal precision with which the BAO scale could be 
measured. 

A detailed analysis of the corresponding constraints on 
the dark energy equation of state i s non-trivial and dep ends 
on the functional form assumed (|Wvithe et al.ll2008h . We 
do not explore the translation of BAO scale to cosmological 
constraints in this paper. Such an analysis should include in- 
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dependent complementary constraints (e.g. from the CMB), 
and take acco unt of uncertainties in, for example, th e hori- 
zon scale (e.g. iMcQuinn et al.ll2006l ; IMao et alj|200gf ). Sim- 
ple estimates of the cosmological constraint s from a high 
redshi ft BAO measurement were presented in IWvithe et al.l 
(|2008h . 

In a recent paper IMao et ah! (I2008T I propose a 
parametrization of the 3d 21 cm PS including peculiar veloc- 
ity effects which is motivated by the generic form expected 
for the ionisation PS and ionisation-matter PS. Our analysis 
of the accuracy with which the BAO scale may be recov- 
ered in the presence of a scale dependent 21 cm bias is com- 
plementary to their detailed discussion of the c osmological 
constr aints available via the redshifted 21 cm PS. lMao et al.l 
(2008) demonstrate that their parametrization can recover 
the ionisation PS constructed from simulations of reionisa- 
tion and therefore can be used to constrain cosmological 
parameters from the 21 cm PS. Whilst they mention that 
their parametrization struggles to describe the high-fc ion- 
isation PS by a redshift of z — 7 when the ionisation PS 
become significantly scale dependent, this is not found to be 
a problem up to the fc max = 0.2 Mpc -1 that we consider in 
this paper. It would be interesting to see how well such a 
parametrization would perform at the end of reionisation in 
the wavelength range that we have used for our fits. 

Our calculations of the constraints on the BAO scale 
from 21 cm intensity fluctuations are promising when com- 
pared to the expecte d performance of future galaxy surveys. 
iGlazebrook fc Blakd l|2005l ) show that a dedicated next- 
generation spectroscopic galaxy survey could measure the 
BAO scale to an accuracy of ~ 1 per cent for the transverse 
component at a redshift z ~ 3.5, with a slightly reduced 
sensitivity to the line-of-sight BAO scale. Limitations in the 
precision achievable arise due to the large number of galax- 
ies needed to reduce the cosmic variance on the scales at 
which the BAO signal is expected (> 30 Mpc), combined 
with the decrease in the number of sufficiently luminous 
galaxies toward higher redshift. The accuracy with which 
the BAO scale may be measured via the 21 cm intensity 
PS is immune to both of these constraints because there 
is no need to resolve individual objects. Redshifted 21cm 
measurements of the PS will instead be limited by cosmic 
variance owing to the finite size of the field of view, and 
sensitivity to individual modes owing to array configuration 
and finite collecting area. However we find that the PS of 
redshifted 21 cm fluctuations as measured by the MWA5000 
could achieve constraints on the BAO scale at z > 6 that 
are comparable to the accuracy which will be available to 
future galaxy redshift surveys at z < 3.5. 

In summary, one potential obstacle to measuring the 
BAO scale in 21 cm fluctuations during the reionisation era 
is the strongly scale dependent bias between the 21 cm and 
mass PS that is produced by the appearance of large ionised 
regions. In this paper we have shown that the appearance 
of ionised regions during the reionisation epoch will not in- 
hibit the measurement of the BAO scale during most of the 
reionisation era. 
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